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We discuss a number of cold nuclear matter effects that can modify 
open heavy flavor and quarkonium production in proton-nucleus collisions 
and could thus also affect their production in nucleus-nucleus collisions, 
in addition to hot quark-gluon plasma production. We show some results 
for p+Pb collisions at yfs = 5 TeV at the LHC. 
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High energy heavy-ion collisions, especially those at collider energies at RHIC 
and the LHC, study the low baryon density and high temperature region of the 
phase diagram of nuclear matter. These experiments are designed to search for phase 
transitions from normal matter to a deconfined quark-gluon plasma. There are other 
mechanisms that can mimic some of the effects of the deconfinement phase transition 
that do not involve quark-gluon plasma production. The way to separate these ’cold 
nuclear matter’ effects from those of the deconfined media is to study other, smaller 
collision systems. Thus a complete study of quark-gluon plasma production and 
evolution includes proton-proton and proton-nucleus collisions at, if not the same 
energy, as close as possible. At RHIC, where the beams are independent, pp, d+Au 
and Au+Au have been studied at the same nucleon-nucleon center of mass energy, 
A ys NN = 200 GeV, as well as other, lower energies. The LHC collider has studied is 
set up so that the beams are not independent, thus the intermediate p+Pb collisions 
run in 2013 involved a 4 TeV proton beam on a 4(Zpb/Apb) = 1.58 TeV lead beam. 
The resulting s NN = 5.02 TeV p+Pb system is close in energy to the upcoming 
5.1 TeV Pb+Pb collisions in LHC Run 2. However, the center of rapidity in these 
collisions is shifted backward by Ay = 0.46 units relative to the symmetric pp and 
Pb+Pb collisions run in Run 1 at yfs = 2.76, 7 and 8 TeV (pp) and = 2.76 

TeV with forward y defined as the direction of the proton beam. In addition, the 
LHCb and ALICE detectors which have muon spectrometers on only one side of the 
detector require two modes of running, one in which the proton beam travels toward 
forward rapidity and one in which it travels toward negative rapidity. 

Deviations from the pp baseline in p+Pb and Pb+Pb collisions have been typically 
been quantified by ratios of observables in the heavy systems to the light system, such 
as the nuclear modification factor, 

D / \ da pPh (y,p T )/dydp T 

RpMy,PT) = Tp —t - 7 - w, T • ( 1 ) 

T pPh dcr pp (y, p T )/dydp T 

So far, no pp collisions have been run at the same energy. Thus the p+Pb data have 
been compared to an interpolation of the pp results. An alternate way of separating 
the cold matter effects in these collisions is the use of the forward-backward ratio, 


da pPh (y > 0 ,p T )/dydp T 
FB V)Pt da pPh (y < 0 ,pr)/dydpr 

which requires no interpolation and cancels systematic ratios in the same system. 

Interpreting the results of these collisions requires an understanding of the cold 
matter effects which are important in their own right. The focus has been on hard 
probes (with a large mass or transverse momentum scale) of these collisions which, 
although produced early in the collision history, can be influenced by the bulk which 
is assumed to reach thermal equilibrium so that its evolution can be described by 
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hydrodynamics. A special advantage of hard probes is that their production, in 
pp collisions, and with suitable modifications, as discussed in the remainder of these 
proceedings, can be described by hard scattering in perturbative QCD. Already nearly 
thirty years ago, J/-0 production was predicted to be suppressed by deconfinement 
in nucleus-nucleus collisions [1]. It soon became clear that suppression of a sort was 
already present in fixed-target proton-nucleus collisions and whatever was causing 
this suppression may also be causing at least some of the suppression in heavy-ion 
collisions. Since the J/ijj is more easily detected by a peak in the dilepton continuum 
at 3.097 GeV, it has been studied longer than open charm hadrons in the same 
systems. The best studied charm mesons are D°jD since they primarily decay into 
two light hadrons, K ± n T , and can thus be reconstructed, even in the high multiplicity 
environment of high-energy heavy-ion collisions. However, the large energies of the 
LHC allow higher statistics measurements of both open and hidden charm, permitting 
detailed studies of other charm mesons and an extension of the charm studies to 
heavier bottom. 

In this brief overview, we will describe some of the most studied cold nuclear mat¬ 
ter effects. These include: isospin (the difference between neutrons and protons as 
projectiles); modifications of the parton distributions in the nuclear medium (shad¬ 
owing); energy loss as the produced cc pair travels through the medium (which can 
manifest itself as a shift in the longitudinal momentum or as a broadening of the 
transverse momentum distribution); intrinsic charm in the proton; and, more spe¬ 
cific to quarkonium final states absorption by nucleons and breakup by interactions 
with produced, comoving particles. We will briefly touch upon all of these except for 
intrinsic charm, covered in a recent review [2]. 

Isospin affects all perturbative QCD calculations because a nucleus is not made up 
of all protons or all neutrons. It is negligible for J/-0 in most regions of phase space 
because J/ip production is dominated by gg processes except in the most forward 
regions. However, it is very important for Drcll-Yan, direct photon and gauge boson 
production. 

The modification of the parton distributions in nuclei, often referred to as shadow¬ 
ing, also affects production of all hard probes. Since nuclear deep inelastic scattering 
experiments probe only charged parton densities, the nature and magnitude of the 
effect on the nuclear gluon density is not well known. The effects of shadowing in nu¬ 
clei are parameterized by various groups using global fitting methods similar to those 
used to evaluate the parton densities in the proton. Currently, LO and NLO sets are 
available, evolving quarks, antiquarks and gluons separately with Q 2 . For example, 
the EPS09 release [3] includes 31 error sets. The EPS09 uncertainty band is obtained 
by calculating the deviations from the central value for the 15 parameter variations 
on either side of the central set and adding them in quadrature. The results for the 
EPS09 LO and NLO gluon modifications at the J/ijj production scale are shown in 
Fig. Ufa). There are considerable differences between the LO and NLO extractions 
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which, however, do not need to manifest themselves as different observables. Indeed, 
in a consistent global analysis, the LO and NLO shadowing ratios for J/ip should 
agree. This is not quite true for the EPS09 set, as shown in Fig. [Q(b) for R p pb{y) 
at 5.02 TeV. The differences between the shadowing ratios calculated for J/ijj pro¬ 
duction with the LO set at leading order and the NLO set at next-to-leading order 
demonstrate this fact. This difference is not an artifact of the more complex NLO 
process since applying the NLO set to the LO calculation gives a result quite similar 
to the consistent NLO calculation. The variation in the same can be attributed to 
the 2—^3 processes dominating cc production at high energies rather than the 2 —» 2 
in the LO color evaporation model. The scales are also somewhat different in the two 
calculations since the NLO calculation in integrated over the full px range. For more 
details, see Ref. [2]. 



Figure 1: (Color online) The EPS09 LO (blue bands) and NLO (red bands) nuclear 
gluon distributions for a lead nucleus, A = 208, are presented as a function of x for 
Q 2 = m 2 (a). Panel (b) shows calculations of the J/-0 nuclear modification factor 
in p+Pb collisions at ^s NN = 5.02 TeV compared to results with EPS09 NLO in a 
NLO calculation (red); EPS09 LO in a fully LO calculation (blue); and EPS09 NLO 
and a fully LO calculation (black) compared to the ALICE [5j !6] and LHCb (?j data. 
(Note that the scale values for m,j are the values for the central renormalization 
and factorization scales from the open charm fits [8].) 

The dependence of nuclear shadowing on collision centrality (central collisions 
have almost full overlap of the two nuclei while more peripheral collisions have smaller 
overlap) has recently been studied in heavy-ion collisions. The PHENIX d+Au J/rj) 
data suggested a stronger dependence than expected from a linear dependence on the 
path length of the probe through the nucleus [9]. We checked whether some power n 
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of the nuclear thickness function, Ta(tt), at impact parameter rp, [10 

= 1 - (1 - R a (x,Q 2 ) ■ 


( 3 ) 


could describe the data. Here R g (x,Q 2 ) is the EPS09 NLO central set and a(n) is a 
normalization adjusted to give the average EPS09 modification when integrated over 
Tt- The PHENIX data preferred n > 10 in Eq. (JUT) . Since this result is physically 
unsatisfying, also We found that assuming shadowing is concentrated in the core 
of the nucleus provides a good description of the data. The formulation for gluon 
shadowing in the core of the nucleus is [ ID] 


K re (r T ) 


1 - 


Rg{x,Q 2 ) 


a(R, d)( 1 + exp((r T — R)/d )) 


(4) 


Here the normalization a(R, d ) is again adjusted to give the appropriate average 
EPS09 modification. The core region is found to be defined by a radius of R ~ 2.4 fm 
with a relatively sharp cutoff, a width of d ~ 0.12 fm. Work on a similar calculation for 
J/ip production as a function of rapidity is in progress for the LHC p+Pb data. More 
data on more final states are needed to test these results. The result for the EPS09s 
spatially-dependent modifications, which retain up to quartic powers in the path 
length of the A-independent coefficients HU, has a much weaker centrality dependence 
than either of the fits to Eqs. (J3]) and (HJ), 

Initial-state energy loss has been studied in Drell-Yan production and found to 
be small. However, there is an inherent ambiguity when applying initial-state energy 
loss to Drell-Yan production since most groups parameterizing the nuclear parton 
densities include these same Drell-Yan data to extract the strength of shadowing on 
the antiquark densities [3]. Also, by forcing the loss to be large enough to explain 
the high xp behavior of J/ip production in fixed-target interactions p2] violates the 
upper bound on energy loss established by small angle forward gluon emission [IB] . 
More recently, it has been proposed that cold matter energy loss should be treated as 
a final-state effect [2]. The final-state J/ip energy loss in pA collisions is currently 
implemented as a probability distribution dependent on the energy loss parameter. 
The effect modifies the xf and pr distributions in a rather crude fashion since the 
quarkonium distribution in pp collisions is parameterized as a convolution of factorized 
power laws, oc (1 — x) n (pl/(pl +Pr)) m i rather than using a quarkonium production 

model mm- 

Initial-state energy loss in the medium can be connected to transverse momentum, 
k t, broadening in nuclei relative to pp collisions, also known as the Cronin effect [16] , 

Final-state nuclear absorption, which affects only quarkonium states, involves 
CNM breakup of the (proto)quarkonium state. Absorption is related to the size 
and production mechanism of the interacting state and can be described by a sur¬ 
vival probability, Sj| bs = exp{— J z °° dz'pA^b, z')a£ hs (z — z')} where z' is the production 
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Figure 2: (Color online) (a) Schematic shift in xp distribution caused by energy loss, 
(b) The LHC J/-0 R p pb(y) data from ALICE |5[[6] and LHCb \J\ compared to energy 
loss model of Arleo et al. US ap¬ 


point and z is the dissociation point; pA(b, z') is the nuclear matter density; and cr^bs 
is the effective absorption cross section for quarkonium state C [T7]. The J/ifj has 
been most studied. Larger effects, at least at midrapidity, have been seen for the 
[T2, IB]- The A dependence of the Xc is still largely unknown. The PHENIX \c result 
[TS] has large uncertainties. 

Previous studies have shown the absorption cross section to depend on rapidity (or 
xf)- Increased effective absorption at backward rapidity may be due to interaction or 
conversion inside the target while increased effective absorption at forward rapidity 
may be due to energy loss. The effective J/^ absorption cross section obtained 
from the PHENIX data, in combination with the EPS09 shadowing parameterization, 
exhibits a somewhat parabolic shape, with a minimum at midrapidity and increasing 
at forward and backward rapidity, according to the intuition above. This result was 
found to be independent of the chosen centrality dependence |10j . This result is 
in agreement with the energy dependence obtained at midrapidity in Ref. |19], the 
best fit decreases with increasing y/s NN . In the same paper, the increase of the 
effective absorption cross section for y > 0, turns on closer to y — 0 for lower y/s NN: 
in addition to the stronger absorption at lower energies | T9] , 

Intrinsic charm is still controversial but could play a role in the dependence of 
forward J/^ production [2DJ (2T|. See Ref. [2] for a recent review and Ref. [22] on the 
interpretation of global fits of the charm structure function. 

So far, only charmonium production has been discussed in detail. Of the afore¬ 
mentioned effects, only absorption by nucleons and comovers is ineffective for open 
charm. The effect of kp broadening was applied to J/i/j production to take the place 
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of low Pt resummation on the cc pair px distribution in pp collisions. No modifica¬ 
tion was assumed in pA collisions. The kx broadening applied to open charm was 
assumed to be the same as in J/ip. In addition, fragmentation must also be taken 
into account. Previously kx broadening and Peterson fragmentation were employed 
to explain the px dependence of D meson production in pA and 7 tA interactions [23]. 
The strong momentum shift of the Peterson parameter, e = 0.06, required a relatively 
large k^ to counteract the fragmentation and bring the result into agreement with 
the data, indeed making the combined effect similar to the original charm quark dis¬ 
tribution itself. More recently, the fragmentation function employed in the FONLL 
approach, without any kx broadening, was shown to be softer El- This behavior can 
be reproduced in the MNR code [25] with a reduced Peterson parameter, e = 0.008. 

At collider energies one might expect that the effect of kx broadening is dimin¬ 
ished. While this is certainly the case at high px, in the low px domain of the ALICE 
measurements, it is still important. To investigate the effects of fragmentation and 
kx broadening on the nuclear modification factor, in Fig. [3] we show the calculated 
p+Pb to pp ratios compared to the ALICE data [26], averaged over D °, D + and D* 
and their antiparticles. The same conditions are used for both p+Pb and pp collisions 
(aside from the EPS09 central set employed in p+Pb) for all histograms except the 
black dotted histogram where a large kx , as might be expected, is used in p+Pb colli¬ 
sions. This result shows that, as long as other effects are kept constant in vacuum and 
in media, shadowing dominates and all the values of R p pb(pT ) are consistent. How¬ 
ever, assuming kx broadening increases in medium counters the shadowing effect in 
the ratio. The large uncertainties in the data are consistent with all the calculations, 
including the one with enhanced broadening. 

In conclusion, numerous cold nuclear matter effects have been postulated. The 
most studied, nuclear shadowing, is important, but its strength is unclear. Fixed- 
target results suggest other effects, including energy loss or intrinsic charm, could be 
more important in the forward region. Clear evidence for intrinsic charm requires 
more data focused in the x and Q 2 region where it should dominate. The effects of 
energy loss in cold matter described here should be verified with other processes. A 
new global analysis of all data, including those of the LHC, should be done to provide 
a new shadowing benchmark. 
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